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Abstract 

The Bohr Hamiltonian describing the collective motion of atomic nuclei is modified by 
allowing the mass to depend on the nuclear deformation. Exact analytical expressions are 
derived for spectra and wave functions in the case of a 7-unstable Davidson potential, using 
techniques of supersymmetric quantum mechanics. Numerical results in the Xe-Ba region 
are discussed. 

1 Introduction 

The Bohr Hamiltonian pQ and its extension, the geometrical collective model [21 E], have 
provided for several decades a sound framework for understanding the collective behaviour 
of atomic nuclei. A puzzle remaining unsolved since the early days of the Bohr Hamiltonian 
regards the behaviour of the moments of inertia of atomic nuclei They are predicted 
to increase proportionally to (3 2 , where (3 is the collective variable corresponding to nuclear 
deformation, while experimentally a much more moderate increase is observed, especially for 
well deformed nuclei. In addition, the use of a constant mass has been recently questioned 
[3j, pointing out that the mass tensor of the collective Hamiltonian cannot be considered 
as a constant and should be taken as a function of the collective coordinates. 

On the other hand, the algebraic framework of the Interacting Boson Model (IBM) [5] 
has been very useful in the systematic study of nuclei corresponding to its three limiting 
symmetries [vibrational U(5), axially symmetric deformed SU(3), 7-unstable 0(6)], as well 
as to intermediate cases. In the geometrical limit of the IBM [6 J , obtained through the use 
of coherent states [6], it is worth remembering that in addition to the usual term of the 
kinetic energy, 7r 2 , terms of the form f3 2 ir 2 appear in the 0(6) limiting symmetry and in the 
U(5)-0(6) transition region, while more complicated terms appear in the SU(3) limiting 
symmetry, as well as in the U(5)-SU(3) and SU(3)-0(6) transition regions [7j. These terms 
indicate that it might be appropriate to search for a modified form of the Bohr Hamiltonian, 



1 



in which the kinetic energy term will be modified by terms containing (3 2 , and even by more 
involved terms. 

The above reasoning leads to the consideration of a Bohr Hamiltonian with a mass 
depending on the collective variable (3. Position-dependent effective masses have been con- 
sidered recently in a general framework [S] , demonstrating the equivalence of this approach 
to the consideration of deformed canonical commutation relations, as well as to the consid- 
eration of curved spaces. Furthermore, several Hamiltonians known to be soluble through 
techniques of supersymmetric quantum mechanics (SUSYQM) [9], [10], have been appro- 
priately generalized [TT] to include position-dependent effective masses, the 3-dimensional 
harmonic oscillator being among them [TTj . 

In the present work we are going to show that the Bohr Hamiltonian with a harmonic 
oscillator potential in f3 (to which a term proportional to 1/ ' (3 2 can be added at no cost) can 
be generalized in order to include a mass depending on f3, B — B /(l + a/3 2 ) 2 , where B 
and a are constants. Exact solutions will be constructed using techniques of SUSYQM. In 
order to achieve exact separation of variables, we are going to limit ourselves to potentials 
independent of the collective variable 7, called 7-unstable potentials. Numerical results for 
the Xe-Ba isotopes, well known examples of 7-unstable behaviour [12], will also be shown. 



2 Formalism of position-dependent effective masses 

When the mass m(x) is position dependent [5], it does not commute with the momentum 
p = —ihV. As a consequence, there are many ways to generalize the usual form of the 
kinetic energy, p 2 /(2mo), where tuq is a constant mass, in order to obtain a Hermitian 
operator. In order to avoid any specific choices, one can use the general two-parameter 
form proposed by von Roos [13], with a Hamiltonian 
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[m d '(x)Vm K '(x)Vm A '(x) +m A '(x)Vm K '(x)Vm d '(x)] + V{x) 
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where V is the relevant potential and the parameters 5', n', X' are constrained by the 
condition 5' + k' + A' = — 1. Assuming a position dependent mass of the form 



m(x) = m M(x), M(x) 



(/(x)) s 



/(x) = l + </(x), 



(2) 



where m is a constant mass and M(x) is a dimensionless position-dependent mass, the 
Hamiltonian becomes 

H = --^[/(x)Vf(x)V/ A (x) + / A (x)Vr(x)V/ 5 (x)] + V{x), (3) 

47Bo 

with 5 + k + X = 2. It is known [S] that this Hamiltonian can be put into the form 
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V7wv/(x)vv / 7(xT+ K//(x), 



with 



V^,(x) = V(x) + 



2m 



[1-6- A)/(x)V 2 /(x) + 



A (V/(x))< 



(4) 
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3 Bohr Hamiltonian with deformation-dependent ef- 
fective mass 



The original Bohr Hamiltonian pQ is 



H 



IB 



1 d 4 d 1 d d 

ft*dft dft + /5 2 sin3 7 ^ Sm 4ft 2 — snr 



--T- 



Q 



(7 - fvrA;) 



(6) 

where /3 and 7 are the usual collective coordinates (ft being a deformation coordinate 
measuring departure from spherical shape, and 7 being an angle measuring departure from 
axial symmetry), while Qk (k = 1, 2, 3) are the components of angular momentum in the 
intrinsic frame, and B is the mass parameter, which is usually considered constant. 

We wish to construct a Bohr equation with a mass depending on the deformation coor- 
dinate ft, in accordance with the formalism described above, 



B(ft) 



Bn 



where Bq is a constant. We then need the usual Pauli-Podolsky prescription [13] 



(7) 



(v$y = g* 



<9$ 



V 2 $ 



in order to construct a Schrodinger equation of the form of Eq. Q in a 5-dimensional 
space equipped with the Bohr- Wheeler coordinates ft, 7. Since the deformation function / 
depends only on the radial coordinate ft, only the ft part of the resulting equation will be 
affected, the final result reading 



2 /3 4 dft P J dft 



f 



P d d 
2ft 2 sin 37 c?7 c>7 



E 



Ql 



8 ^,3 sin2 (7-§^) 



+ V eff 



^ = e^, (9) 



where reduced energies e = B E/h 2 and reduced potentials v = B V/h 2 have been used, 
with 

v eff = v(ft tl ) + -(1 - 5 - A)/V 2 / + 1 / ' A ' ' 
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A (V/)' 
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4 The 7-unstable Davidson potential 

The solution of the above Bohr-like equation can be reached for certain classes of potentials 
using techniques developed in the context of SUSYQM [9j [101 [11]. In order to achieve 
separation of variables we assume that the potential v(ft, r j) depends only on the variable 
ft, i.e. v(ft,j) = u(ft) [15J. Potentials of this kind are called 7-unstable potentials, since 
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they are appropriate for the description of nuclei which can depart from axial symmetry 
without any energy cost. Furthermore, we are going to use as an example the Davidson 
potential [TO] 



u{(3) =(3 2 + 



ll 2 
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where the parameter /3o indicates the position of the minimum of the potential, the special 
case of /3q — corresponding to the simple harmonic oscillator. (Note that the term 
containing /3 offers no additional complication to the solution). 
One then seeks wave functions of the form [T~5l ITT] 



where 6; 



m3,-y,0 i ) = F(/3)$(y,e i ), 
1, 2, 3) are the Euler angles. Separation of variables gives 



(12) 



2 I3 A d(3 P 1 d(3 
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A)/V 2 / 

- X] (V/) 2 + u({3) 



F(P) = eF(fi), (13) 
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(7 - Ink) 
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$( 7 ,^)=A$( 7 ,^ 



(14) 



Eq. (14) has been solved by Bes [18] . A = r(r + 3) represents the eigenvalues of the second 
order Casimir operator of SO (5), while r is the seniority quantum number, characterizing 
the irreducible representations of SO (5). The values of angular momentum L occurring 
for each r are provided by a well known algorithm and are listed in [HI EES]. Within the 
ground state band (gsb) one has L = 2r. The L = 2 member of the quasi-71 band is 
degenerate with the L = 4 member of the gsb, the L = 3, 4 members of the quasi-71 band 
are degenerate to the L = 6 member of the gsb, the L = 5, 6 members of the quasi-71 band 
are degenerate to the L = 8 member of the gsb, and so on. 
Eq. (13) can be simplified by performing the derivations 



-fF" + iff + 



2f 



F' + 



if) 2 | // 
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P 
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F - - — AF + eF 

2(3 2 



v eff F = 0, (15) 



with 



s 



A {ff- 



(16) 



The difference in the numerical coefficient of /' observed in comparison to Eq. (2.27) of 
Ref. [8] is due to the different dimensionality of the space used in each case. 
Setting 



Eq. ( 15 ) is put into the form 

HR = 



ff 2 



fjpVf) R + u eff R = eR 



(17) 



(18) 
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where 



U eff = Veff + ~ + ^A. (19) 

Based on the results for the 3-dimensional harmonic oscillator reported in Ref. [TT], in 
order to find analytical results for Eq. (18) we are going to consider for the deformation 
function the special form 

f((3) = l + a(3 2 . (20) 

This choice is made in order to lead to an exact solution. Its physical implications will be 
discussed in Section 8. 

Using these forms for the potential and the deformation function one obtains 



u 



eff 



l + a 2 /? 2 



§(l-*-A) + 2(i-* 



1 \ A 
--A +3 + - 

2 J 2 



+ 



1 

2' 



p l l + ^A + /5 4 ) +a 



;i-5-A)+4 + A 



• (21) 



5 Factorization 



Following the general method used in supersymmetric quantum mechanics (SUSYQM) 
(9J [TO] , one should take the following steps: 

i) Factorize the Hamiltonian. 

ii) Write a hierarchy of Hamiltonians starting from the first one. 

iii) Impose the shape invariance conditions, which are integrability conditions guaran- 
teeing exact solvability. 



Thus one first tries to put the Hamiltonian in Eq. ( 18 ) in the form 

Hq = BqBq + e . 



(22) 



Then this Hamiltonian can be considered as the first member of a hierarchy of Hamiltonians 

i 

H % = BlB; + Y;Z3, (23) 

3=0 

expressed in terms of the generalized ladder operators , B^~, which are determined 
recursively, along with e$ [which are in general different from the energy eigenvalues e 
appearing in Eq. (18), the only exception being e = £o], through the shape invariance (SI) 
condition 

B7B+ = B+ +1 B- +1 + e i+l . (24) 
In the present case one can start with 



(25) 
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where the second term resembles the superpotential occurring in the case of the 3-dimensional 

;p< 
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harmonic oscillator, reported in Ref. [IT] . Substituting into Eq. (22) one obtains 
-"0 



dp 



+ 



2[3 2 



^c + ^ac + c c + £ ) ■ ( -0 ) 



Comparing this result to Eq. (|18|) with the effective potential of Eq. (21) and equating 

1 



powers of (3 one gets the following: 

Co 



± y/a? + 8Pi) 



1 ± a/9 + 4A + 8^), and 



so 
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-c - -ac - c c + -(1 - 5 - X)a + Aa + aA, 



where 



P x = 1 + a 2 



A +3 + 



A 



(27) 

(28) 
(29) 

(30) 



Upon substitution of c and c from Eqs. (27) and (28), Eq. (29) leads to 
'29 



£o 



4 



--(5 + X)+A ± 



1 



- x /9 + 4A + 8/3 4 



- C r- x /(a 2 + 8Pi)(9 + 4A + 8/3^ 



(31) 



where ± refers to the sign in Eq. (27), =F corresponds to the sign selection in Eq. (28), 
while a is +1 for the selections (+, — ) and (— , +) in Eqs. (27) and (28) respectively, while 
it is —1 for the selections (+, +) and (— , — ). 

The limiting case of a = 0, in which Pi — 1, can be used for reducing the number 
of possibilities. Since Eq should be increasing with A (members of the ground state band 
should exhibit energies increasing with the angular momentum L), it turns out that one 
should have a = +1, thus only the choices (+, — ) and (— , +) are acceptable. 

In Section [7] it will be shown that good behaviour of the wavefunctions at /3 = leaves 
(+, — ) as the only choice. 



6 Shape invar iance 



In the next step the hierarchy of Hamiltonians of Eq. ( 23 ) should be considered, with 

d 



Bf = T 



1 

71 



d/3 



f + 



V2\ 







(32) 



Substituting these expressions in the shape invariance condition of Eq. (24) and equating 
powers of (3 leads to the following results: 



2 , _ 2 

Cj -+- CjCl — Cj +1 — Q+id, 



(33) 
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$ ~Ci = c-+i + Q+i, (34) 

and 

2£j+i = Q + Cj + i — aQ — acj + i + 2QCj — 2cj + iCj + i. (35) 

Keeping from the first two of these oniy the solutions c i+ \ = Ci + a (leading to q = c + ia) 
and Cj+i = Cj — 1 (leading to Cj = Cq — i), in accordance with the results obtained for the 
3- dimensional harmonic oscillator [11], we get 

Si = 2[c — ac + a(2i — 1)]. (36) 

One then easily finds for the energy 

e„ = Si — e + 2z/(c - ac ) + lav 2 . (37) 

The ground state band is obtained for v — 0, while the quasi- /5i band corresponds to v = 1. 



Equation (37) only provides a formal solution to the bound-state energy spectrum. 
The range of v values is actually determined by the existence of corresponding physically 
acceptable wavef unctions, to be discussed in the next Section. 

7 Wave functions 

To be physically acceptable, the bound-state wavefunctions should satisfy two conditions 
P: 

(i) As in conventional (constant-mass) quantum mechanics, they should be square integrable 
on the interval of definition of u e g, i.e., 

df3\R 1/ (/3)\ 2 < oo. (38) 

(ii) Furthermore, they should ensure the Hermiticity of H. For such a purpose, it is enough 
to impose that the operator ^ff{d/df3)\ff be Hermitian, which amounts to the restriction 

\R»(P)\ 2 f([3) -> for (3 -f and (3 -> oo, (39) 

or, equivalently, 

\R v (/3)\ 2 ^0 for/3->0 and \R U {(3)\ 2 (3 2 -> for /3 -> oo. (40) 



As condition (40) is more stringent than condition (38), we should only be concerned with 
the former. 



The ground state wave function can be determined from the differential equation [TT] 

B (p;c ,Co)Ro(P;c ,c ) = 0, (41) 



where B is given by Eq. (25). Trying the solution 

#o = C p n f, (42) 
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where Co is a normalization constant, the powers of (3 lead to the conditions 

n = -c , 

and 



n 



aco — cq — a 



0, the function \R 



2a 

2 behaves as (3 2n . 



(43) 



(44) 



Condition (40) imposes that n > 0, 



For (3 

i.e. Cq < 0. From Eq. (|28l) it is clear that this is guaranteed if the minus sign is retained 



in it. From the discussion given at the end of Section [5] it is clear that the only possibility 
remaining for the signs in Eqs. (27) and (28) is the (+, — ) one. One can easily see that 
this choice guarantees n < 0. 

For (3 — > oo, \Ro{/3)\ 2 f3 2 behaves as (3~ 2c °^ a , since in this case / rs a/3 2 . Condition (40) 



therefore imposes that Co > 0. This restriction is already satisfied, since we have been led 
in the previous paragraph to keep the upper sign choice in (27). 



Wave functions of excited states can then be obtained from the recursion relation [TTj 

1 



R„ +1 (f3; c , Co) 



e u+ i(co,c ) - e (co,c J 



-.B + (/3]Co 1 co)R„(/3;c 1 ,c 1 ), 



where the different coefficients appearing in the last term should be noticed, 
recursion relation one obtains 



(45) 



From the 



Ri{(3] c , c ) 



Ci 



(2c + a)f3 + 



2c n 



1 







n+l rn—1 



P U+ f 



(46) 



# 2 (/3;c ,c ) 



Co 



2 3 /V(e 2 -e )(ei-e ) 



(2c + 3a)(2c + a)l3 2 + (2c - 3)(2c - + 2(2c + 3a)(2c - 3) 



n+2 ffi-2 



P n+2 f 



(47) 



where C\, Ci are normalization constants. 

From the above it is clear that wave functions of the states belonging to the ground state 
band are obtained from Eq. (12), substituting in it Eqs. (17) and (42), while wave functions 



of the states belonging to the quasi- {3% band are obtained from Eq. (12), substituting in it 
Eqs. fll7} and dil. 



It is easy to see that the conditions imposed above in order to guarantee the physi- 
cally acceptable behaviour of the ground state wavef unctions, also guarantee the physically 
acceptable behaviour of the wavefunctions for excited states. In particular: 

(i) In order to examine the behaviour at (3 — > 0, it suffices to examine in R v the behaviour 
of the polynomial term containing the lowest power of (3 {f3~ l in Ri, (3~ 2 in R 2 ). In both 
cases the function \R v (f3)\ 2 behaves as /3 2n , i.e. it exhibits the same behaviour as \R (j3)\ 2 . 

(ii) In order to examine the behaviour at (3 — > oo, it suffices to examine in R v the behaviour 
of the polynomial term containing the highest power of (3 ((3 in Ri, f3 2 in R 2 ). In both cases 
the function \R V {(3)\ 2 (3 2 behaves as /3~ 2c0//a , i.e. it exhibits the same behaviour as \Ro{(3)\ 2 (3 2 . 

Therefore the wavefunctions of the excited states given above are forced to exhibit 
physically acceptable behaviour by the same conditions which guarantee the physically 
acceptable behaviour of the ground state wavefunctions. 
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8 Numerical results 



FromEq. Q it is clear that in the present case the moments of inertia are not proportional 
to (3 2 sin 2 (7 - 27rk/3) but to {(3 2 /f 2 {(3)) sin 2 (7 - 2vr£;/3). The function (3 2 /f 2 {(3) is shown 
in Fig. 1 for different values of the parameter a. It is clear that the increase of the moment 
of inertia is slowed down by the function f(/3), as it is expected as nuclear deformation sets 
in 0| . 

As a first testground of the present method we have used the Xe isotopes shown in Table 
1. They have been chosen because: 

i) They are known to lie in a 7-unstable region [T2] . 

ii) At least the bandheads of the quasi-/?i and quasi-71 bands are known experimentally. 

iii) They extend from the borders of the neutron shell ( 134 Xeso is just below the N=82 
shell closure) to the midshell ( 120 Xe66) and even beyond, exhibiting increasing collectivity 
(increasing = E(4f)/E(2f) ratios) from the border to the mishell. 

For evaluating the rms fits performed, the quality measure 



" = ' ! (n - TO (48) 
has been used. The theoretical predictions for the levels of the ground state band are 



obtained from Eq. (31 ) (in which all terms with double signs are taken with positive signs, 



as explained at the end of Section 5), while the levels of the quasi- (3\ band are obtained from 



Eq. (37) for v — 1. The levels of the quasi-71 band are obtained through their degeneracies 



to members of the ground state band, mentioned below Eq. (14). 

Moving from the border of the neutron shell to the midshell, the following remarks apply 

i) 134 Xe and 132 Xe are almost pure vibrators. Therefore no need for deformation depen- 
dence of the mass exists, the least square fitting leading to a = 0. Furthermore, no /3 term 
is needed in the potential, the fitting therefore leading to fio = 0, i.e., to pure harmonic 
behaviour. 

ii) In the next two isotopes ( 130 Xe and 128 Xe) the need to depart from the pure harmonic 
oscillator becomes clear, the fitting leading therefore to nonzero (3$ values. However, there 
is still no need of dependence of the mass on the deformation, the fitting still leading to 
a = 0. 

iii) Beyond 126 Xe both the (3q term in the potential and the deformation dependence of 
the mass become necessary, leading to nonzero values of both flo and a. 

Exactly the same behaviour is seen in the Ba isotopes, also known to lie in a 7-unstable 
region [12] and shown in Table 1. 

The results shown in Table 1 have been obtained for 5 = A = 0. One can easily verify 
that different choices for 8 and A lead to a renormalization of the parameter values a and 
/3q, the predicted energy levels remaining practically the same. 

In addition to energy spectra, B(E2) transition rates should be calculated and compared 
to experiment. The details of this task are deferred to a longer publication. However, basic 
qualitative features can be seen in Fig. 2, where the systematic behaviour of energy ratios 
and B(E2) ratios within the ground state band are shown. We remark that an increase of 
the a parameter leads to a more rapid increase of energies (normalized to the energy of the 
first excited state) within the ground state band as a function of the angular momentum 
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L, while in parallel it slows down the increase of B(E2)s (normalized to the transition from 
the first excited state to the ground state) as a function of L. 

Among the nuclei of Table 1, the only one exhibiting experimentally known increasing 
B(E2)s within the ground state band is 128 Xe. The theoretical predictions (using the param- 
eters of Table 1, obtained by fitting the energy levels alone) fall withing the experimental 
error bars, as seen in Table 2. 

The present results suggest that dependence of the mass on deformation becomes nec- 
essary as deformation increases. It is therefore desirable to provide a similar solution of 
the Bohr Hamiltonian applicable to axially symmetric well deformed nuclei. Work in this 
direction is in progress. 

9 Conclusion 

Motivated by the existence in the geometrical limit of the 0(6) limiting symmetry and of the 
U(5)-0(6) transition region of the Interacting Boson Model of extra terms of the form /3 2 7r 2 
(where (3 is the nuclear deformation) in addition to the kinetic energy term it 2 , as well as by 
the existence of more complicated additional terms in the SU(3) limiting symmetry and in 
the U(5)-SU(3) and SU(3)-0(6) transition regions, we have modified the Bohr Hamiltonian 
describing the collective motion of atomic nuclei by allowing the mass to depend on the 
nuclear deformation. Using techniques of supersymmetric quantum mechanics we have 
obtained exact analytical expressions for spectra and wave functions for the case of the 
7-unstable Davidson potential. A first numerical application in the Xe-Ba region gives 
encouraging results. Detailed comparisons to experiment, including B(E2) transition rates, 
as well as extension of the method to deformed axial nuclei (with 7^0) are deferred to a 
longer publication. 
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Figure 1: The function /3 2 /f 2 (/3) = (3 2 / (1 + a/3 2 ) 2 , to which moments of inertia are propor- 
tional as seen from Eq. ([9]), plotted as a function of the nuclear deformation (3 for different 
values of the parameter a. See Section 8 for further discussion. 
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Table 1: Comparison of theoretical predictions of the 7- unstable Bohr Hamiltonian with 
/^-dependent mass (with 5 = A = 0) to experimental data [19] of Xe and Ba isotopes. The 
R4/2 = E(4f)/E(2±) ratios, as well as the quasi- (3\ and quasi-71 bandheads, normalized to 
the 2+ state and labelled by R Q/2 = £(0+)/£(2+) and R 2/2 = £(2+)/£(2+) respectively, 
are shown. The angular momenta of the highest levels of the ground state, quasi-/?i and 
quasi-71 bands included in the rms fit are labelled by L g , Lp, and L 7 respectively, while 
n indicates the total number of levels involved in the fit and a is the quality measure of 
Eq. (48). The theoretical predictions are obtained from the formulae mentioned below Eq. 
(48). See Section 8 for further discussion. 
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16 


0.638 


124 Xe 


2.48 


2.36 


3.6 


3.5 


2.4 


2.4 


1.55 


0.051 


20 


2 


11 


21 


0.554 


126 Xe 


2.42 


2.33 


3.4 


3.1 


2.3 


2.3 


1.42 


0.064 


12 


4 


9 


16 


0.584 


128 Xe 


2.33 


2.27 


3.6 


3.5 


2.2 


2.3 


1.42 


0.000 


10 


2 


7 


12 


0.431 


130 Xe 


2.25 


2.21 


3.3 


3.1 


2.1 


2.2 


1.27 


0.000 


14 





5 


11 


0.347 


132 Xe 


2.16 


2.00 


2.8 


2.0 


1.9 


2.0 


0.00 


0.000 


6 





5 


7 


0.467 


134 Xe 


2.04 


2.00 


1.9 


2.0 


1.9 


2.0 


0.00 


0.000 


6 





5 


7 


0.685 


130 Ba 


2.52 


2.42 


3.3 


3.2 


2.5 


2.4 


1.60 


0.118 


12 





6 


11 


0.352 


132 Ba 


2.43 


2.29 


3.2 


2.8 


2.2 


2.3 


1.29 


0.059 


14 





8 


14 


0.619 


134 Ba 


2.32 


2.16 


2.9 


2.7 


1.9 


2.2 


1.12 


0.000 


8 





4 


7 


0.332 


136 Ba 


2.28 


2.00 


1.9 


2.0 


1.9 


2.0 


0.00 


0.000 


6 





2 


4 


0.250 



Table 2: Experimental B(E2;L — > L — 2) transition rates [normalized to the transition 
from the first excited state to the ground state, B(E2;2 — > 0)] within the ground state 
band of 128 Xe [19J, compared to theoretical predictions using the parameters of Table 1. 
See Section 8 for further discussion. 



L 


exp. 


th. 


4 


1.468 ±0.201 


1.632 


6 


1.940 ±0.275 


2.196 


8 


2.388 ±0.398 


2.751 


10 


2.736 ± 1.138 


3.310 
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Figure 2: Energy levels E(L) [normalized to the energy of the first excited state, E(2)] and 
B(E2; L — > L — 2) transition rates [normalized to the transition from the first excited state 
to the ground state, B(E2; 2 — > 0)] are shown for the ground state band as functions of the 
angular momentum L for (3 — 1 and varying values of a (0.0, 0.1, 0.5). See Section 8 for 
further discussion. 
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